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Poisson distribution parameter by single observation of number of events 
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The empirical probability density function for the conditional distribution of the true value of Poisson dis- 
tribution parameter on one measurement is constructed by computer experiment. The analysis of the obtained 
distributions confirms that these distributions are gamma-distributions. 



1. INTRODUCTION 

Let us consider the Gamma-distribution with 
probability density 
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As it follows from the article pQ (see also [2]) 
and is clearly seen from the analysis of identity 3 
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i.e. 
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At change of standard designations of Gamma- 
distribution — , a and x for a, n + 1 and A we 
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get the following formula for probability density 
of Gamma-distribution 
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where a is a scale parameter and n + 1 > is a 
shape parameter. Suppose a = 1, then the proba- 
bility density of Gamma-distribution Pi.„+i looks 
like Poisson distribution of probabilities: 



A™ 

.9«(A) = -re-\ A > 0, n > -1. 
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Let the probability of observing n events in the 
experiment be described by a Poisson distribution 
with parameter A, i.e. 
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for any Ai > and A2 > 0, that at one measure- 
ment of the number of events n (in our case it 
is the number of casual events appearing in some 
system for certain period of time) which appear 
according to Poisson distribution, the parameter 
value of this distribution is described by Gamma- 
distribution r^i+n with mean, mode, and vari- 
ance n + 1, n, and n + 1, respectively. In other 
words conditional distribution of the probability 
of true value of parameter of Poisson distribution 
is a Gamma-distribution T\^ +n on condition that 
the measured value of the number of events is 
equal to n. 

It means that the value n corresponds to the 
most probable value of parameter of Poisson dis- 
tribution, and the mean value of the number of 
events, appearing in Poisson flow in the fixed time 
interval, must correspond to the magnitude n+1, 
i.e. the estimation of parameter of Poisson distri- 
bution at one observation is displaced for 1 from 
the measured value of the number of events. The 
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equation (5) in the considered case allows to mix 
Bayesian and frequentist probabilities. 

As a result, we can easily construct the confi- 
dence intervals, take into account systematics and 
statistical uncertainties of measurements at sta- 
tistical conclusions about the quality of planned 
experiments, estimate the value of the param- 
eter of Poisson distribution by several observa- 
tions 02]. 

Nevertheless there are works in which the ap- 
proaches based on other assumptions of distribu- 
tion of true value of parameter of Poisson distri- 
bution in presence of its estimation on a single 
measurement, for example 5 ]. Also the works 
using Monte Carlo methods for construction of 
confidence intervals and for estimations of Type I 
and Type II errors in the hypotheses testing have 
recently appeared |o17| (see, also, [B]). Therefore 
the experimental test with the purpose to con- 
firm, that the true value of parameter of Poisson 
distribution at single observation has density of 
probability of Gamma-distribution, and with the 
purpose to check up applicability of Monte Carlo 
methods to such tasks was carried out. 

The structure of the paper is the following. In 
the next section the arrangement of computer ex- 
periment is considered, in the third section - the 
statistical analysis of the results is given, and the 
last section contains concluding remarks. 

2. The arrangement of measurements 

From the identity (5) follows, that any prior, 
except the uniform, on value of parameter of Pois- 
son distribution in distribution of true value of 
this parameter at presence of the measured esti- 
mation n is excluded by existence of the boundary 
conditions determined by the appropriate sums 
of Poisson distributions (see Eq.(5)). Therefore 
we carried out the uniform scanning in parameter 
of Poisson distribution with step 0.1 from value 
A = 0.1 up to value A = 20, playing the Poisson 
distribution 30000 trials for each value A (Fig.l) 
with the using of function RNPSSN gj. 

After scanning for each value of number of the 
dropped out events n the empirical density of 
probability of true value of parameter of Poisson 
distribution to be A if the observation is equal n 




Figure 1. Amount of occurrences of n in the interval 
from up to 9. Scanning in parameter of the Poisson 
distribution (30000 trials at each value of parameter 
A) was carried out with step 0.1 in the interval of A 
from 0.1 up to 20 



was obtained. 

3. The analysis of results 

In Fig. 2 the distribution (a), obtained at scan- 
ning (the playing of Poisson distribution with 
consequentive increase of the parameter A after 
each series of trials with the fixed value of the 
parameter) in parameter A with the selection of 
number of the dropped out events n — 6, and 
distribution (b) of the casual value, having Tij 
distribution of appropriate area, calculated by the 
formula, are shown. One can see that the average 
value of parameter A ~ 7. It means, the number 
of observed events is displaced for one unity by 
the estimation of the mean value of Poisson dis- 
tribution parameter and correspond to the most 
probable value (the mean value has bias). 

The same distributions obtained by the selec- 
tion of number of dropped out events n — 
(Fig. 3) and n = 8 (Fig. 4) superimposed on each 
other in logarithmic scale also are shown. 

In Tab.l the values of probabilities of compat- 
ibility of the empirical distribution, obtained by 
Monte Carlo by the scanning in parameter A, and 
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Figure 2. Distributions of occurrences of value n — 6 
depending on value of parameter A. The distribution 
(a) is obtained at Monte Carlo scanning in param- 
eter A. The distribution (b) is obtained by direct 
construction of Gamma-distribution Tij 




Figure 3. Distributions of occurrences of value n — 
depending on value of parameter A. The distribution 
obtained at Monte Carlo scanning in parameter A is 
superimposed on the distribution obtained by direct 
construction of Gamma-distribution Vii 




Figure 4. Distributions of occurrences of value n — 8 
depending on value of parameter A. The distribution 
obtained at Monte Carlo scanning in parameter A is 
superimposed on the distribution obtained by direct 
construction of Gamma-distribution Yi,g 



the appropriate Gamma-distribution for values n 
from up to 9 are presented. The calculations 
are based on the Kolmogorov Test (the function 
HDIFF of the package HBOOK 0). The authors 
of a package as criterion of coincidence of two dis- 
tributions recommend to use the requirement of 
probability value of compatibility more than 0.05. 
In Fig. 4 the least conterminous distributions are 
given. 



Table 1 



The probability of compatibility 



n 


probability 


ii 


probability 





1.000000 


5 


0.999084 


1 


0.999646 


6 


0.999986 


2 


0.992521 


7 


0.999892 


3 


0.999986 


8 


0.752075 


4 


0.999969 


9 


0.974236 



Thus, the obtained results do not contradict 
the statement that conditional distribution of 
true value of parameter of Poisson distribution 
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at single measurement is obeyed to a Gamma- 
distribution. 

4. Conclusion 

In the report Monte Carlo experiment on the 
check of the statement, that true value of parame- 
ter of Poisson distribution at an estimation of this 
parameter on one observation n has probability 
density of Gamma-distribution ri.„+i, is carried 
out. The obtained results confirm the conclusions 
of the papers |3I4| about a kind of conditional 
distribution of true value of parameter of Poisson 
distribution at single observation. 

Note, that the given results also specify the ap- 
plicability of Monte Carlo method for construc- 
tion of conditional distribution of the true value 
of parameters of various distributions. 
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Appendix 

Let us consider the famous equality |1I2| in 
form as written in |3] 

i=n+l "^'1 i=0 

where /xi > 0, /12 > 0, 

m) = f (m; fc ) = -Tj 1 fe = 1, 2, . . . . 

Let us suppose that g(/i; k) is the probability 
density of parameter of the Poisson distribution 
to be \i if number of observed events is equal to k. 
It is a conditional probability density. As shown 
above (Eq.3) the g(n; k) is the density of Gamma- 
distribution by definition. 

On other side: if k) is not equal to this 
probability density and the true probability den- 
sity of the Poisson parameter is the other function 
h(fx; k) (continuous or with set of points of dis- 
continuity by measure 0) then there takes place 
another identity 



00 „ M2 n 

h{fJ-;n)dfi+^2f(i;fi 2 ) = 1, (7) 

i=n+l 111 i=0 

This identity is correct for any [i\ > and 
,x 2 > 0. 

If we subtract Eq.7 from Eq.6 then we have 

/ (g(v, n) - h(p; n))dfj, = 0. (8) 

We can choose the fi% and /i2 by the arbitrary 
way. Let us make this choice so that g(/j,;n) is 
not equal h(n;n) in the interval (/xi,/i 2 ) and, for 
example, g(fj,;n) > h(/J,;n) and ^2 > Mi- in this 
case we have 

I (g(n;n)-h(ji;n))dn>0 (9) 
J in 

and as a result we have contradiction. The iden- 
tity (6) does not leave a place for any prior except 
uniform. 



